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Global perturbative QCD analyses, based on large data sets from electron-proton and hadron col- 
lider experiments, provide tight constraints on the parton distribution function (PDF) in the proton. 
The extension of these analyses to nuclear parton distributions (nPDF) has attracted much interest 
in recent years. nPDFs are needed as benchmarks for the characterization of hot QCD matter in 
nucleus-nucleus collisions, and attract further interest since they may show novel signatures of non- 
linear density-dependent QCD evolution. However, it is not known from first principles whether the 
factorization of long-range phenomena into process-independent parton distribution, which under- 
lies global PDF extractions for the proton, extends to nuclear effects. As a consequence, assessing 
the reliability of nPDFs for benchmark calculations goes beyond testing the numerical accuracy of 
their extraction and requires phenomenological tests of the factorization assumption. Here we argue 
that a proton-nucleus collision program at the LHC would provide a set of measurements allowing 
for unprecedented tests of the factorization assumption underlying global nPDF fits. 



Parton distribution functions (PDFs) fi/h( x ,Q 2 ) play 
a central role in the study of high energy collisions in- 
volving hadronic projectiles h. They define the flux of 
quarks and gluons (i — q, g) in hadrons as a function of 
the partonic resolution scale Q 2 and hadronic momen- 
tum fraction x. For protons, sets of collinearly factorized 
universal PDFs have been obtained, since a long time 
ago, in global perturbative QCD analyses. These are 
based on data from deep-inelastic lepton-proton scatter- 
ing (DIS) and Drell-Yan (DY) production, as well as W/Z 
and jet production at hadron colliders. These data pro- 
vide tight constraints on PDFs over logarithmically wide 
ranges in Q 2 and x, and have allowed precision testing 
of linear perturbative QCD evolution. In comparison to 
proton PDFs, our understanding of parton distribution 
functions fi/A{x,Q 2 ) in nuclei of nucleon number A is 
much less mature. Knowledge of nuclear parton distri- 
bution functions (nPDFs) is important in heavy ion col- 
lisions at RHIC and at the LHC for a quantitative con- 
trol of hard processes, which are employed as probes of 
dense QCD matter. Characterizing nuclear modifications 
of PDFs is also of great interest in its own right, since 
the nuclear environment is expected to enhance parton 
density-dependent effects, which can reveal qualitatively 
novel, non-linear features of QCD evolution. 

Paralleling the determination of proton PDFs, several 
global QCD analyses of nPDFs have been made within 
the last decade [lHal- Up until recently, these analyses 
were based solely on fixed-target nuclear DIS and DY 
data. Compared to the data constraining proton PDFs, 
these are of lower precision and lie in a much more limited 
range of Q 2 and x. Constraints on nuclear gluon distri- 
bution functions are particularly poor, since they cannot 
be obtained from the absolute values of DIS structure 
functions, but only from their logarithmic Q 2 -evolution, 



for which a wide Q 2 -range is mandatory. To improve 
on this deficiency, recent global nPDF analysis [l], Q 
have included for the first time data from inclusive high- 
er hadron production in hadron-nucleus scattering mea- 
sured at RHIC 

However, in contrast to the theoretical basis for global 
analyses of proton PDFs, the separability of nuclear 
effects into process-independent nPDFs and process- 
dependent but A-independent hard processes is not es- 
tablished within the framework of collinear factorized 
QCD. In particular, some of the characteristic nuclear 
dependencies in hadron-nucleus collisions, such as the 
Cronin effect Q, may have a dynamical origin that 
cannot, or can only partly, be absorbed in process- 
independent nPDFs. In view of the importance of nPDFs 
for characterizing benchmark processes in heavy ion colli- 
sions, it is thus desirable to look for stringent phenomeno- 
logical tests of the working assumption of global nPDF 
fits that the dominant nuclear effects can be factorized 
into the incoming PDFs. Here, we argue that a program 
of hadron-nucleus collisions at the LHC would provide 
for such tests with unprecedented quality. 

We will focus mainly on single inclusive high-pT hadron 
production. In the factorized QCD ansatz to hadron- 
nucleus collisions, the cross section for production of a 
hadron h takes the form 

d 3 aP A^h x =AJ2 fj/p fi/A ® d 3 a^ k x ® D k ^ h , (1) 

ijk 

where the symbol £g> stands for the convolution of the in- 
coming PDFs with the cross section of the hard partonic 
process and with the fragmentation function for a parton 
k into a hadron h. The sum goes over all parton species 
contributing to the production of h. By construction, the 
entire nuclear dependence of the cross section (JlJ resides 
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in the nPDF fi/A(x, Q 2 )- It is customary to characterize 
nuclear effects by the ratios 



P T (GeV) 



Rf (x,Q 2 ) = f l/A (x,Q 2 )/f l/p (x,Q 2 



(2) 



In global nPDF analyses, characteristic deviations of 
Rf(x, Q 2 ) from unity are found for all scales of Q 2 tested 
so far and for essentially all scales of the momentum frac- 
tion x. These effects are typically referred to as nuclear 
shadowing (x < 0.01), anti-shadowing (0.01 < x < 0.2), 



EMC effect (0.2 



< 



< 0.7) and Fermi motion (x > 



0.7). A typical example for the nuclear cc-dependence of 
Rf(x, Q 2 ) is shown in the upper left plot of Fig. [TJ 

Nuclear effects on single inclusive hadron production 
are typically characterized by the nuclear modification 
factor RpA, which depends on the transverse momentum 
Pt and the rapidity y of the hadron, 



RpA(PT,y) 



d(J pA^h+X 

dp\ dy 



N. 



pA 



d(jPp ^h+X 

dp\ dy 



(3) 



Here, N^ on denotes the average number of equivalent 
nucleon-nucleon collisions in a pA collision. It is deter- 
mined by Glauber theory, which can be subjected to in- 
dependent phenomenological tests. The lower left plot of 
Fig. [T] shows the nuclear modification factor RjAuiPT, V) 
for the production of neutral pions in ^/snn = 200 
GeV deuteron-gold collisions at RHIC, calculated within 
the factorized ansatz ([1]) at leading order (LO). Re- 
sults shown in Fig. [T] are also consistent with the NLO- 
calculation of Rd°A u (PT, y) m Q- All our calculations use 
LO PDFs from CTEQ6L [l(| with nuclear modifications 
EPS09LO [1] and the KKP fragmentation functions [H|. 
We have checked our conclusions for another set of frag- 
mentation functions [ijj (data not shown). 

The p^-dependence of the nuclear modification factor 
traces the x-dependence of nPDFs. The precise kine- 
matic connection between the momentum fractions xi, 
X2 and the measured hadronic momentum px is com- 
plicated by the convolution of the distributions in (fTJ). 
Qualitatively, at fixed rapidity y of the produced hadron, 
increasing pt tests larger values of x\ , X2 ■ Inspection of 
the nuclear modification factor in the lower left panel 
of Fig. Q] reveals that the enhancement of RdAu(PT,y) 
in the region around pr — 4 GeV at mid-rapidity tests 
momentum fractions in the anti-shadowing region. The 
RHIC data @ in Fig.[T]have been used in constraining the 
nPDF analysis EPS09 [l| but they were not employed in 
a closely related nPDF fit which provides an equally 
satisfactory description of these RHIC data. Therefore, 
the agreement of data and calculation in Fig. Q] is in sup- 
port of collinear factorization. 

However, qualitatively different explanations of the 
RdAu(PT,y) measured at RHIC are conceivable. The 
above calculation accounted for RdAu(PT, V — 0) m terms 
of a nuclear modification of the longitudinal parton mo- 
mentum distribution, only. Alternatively, it has been 
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FIG. 1: (left,top) The ratio © of nuclear to nucleoli PDFS for 
valence up-quarks at Q 2 = (10 GeV) 2 obtained in the EPS09 
LO analysis. Dashed lines characterize the range of uncer- 
tainties, (left, bottom) The nuclear modification factor © for 
neutral pion production in ^/snn = 200 GeV dAu (RHIC). 
Data from PHENIX Q are compared to a EPS09 LO cal- 
culation. Hereafter, uncertainty bands are from EPS09 LO 
only (uncertainties from proton PDFs and FF are neglected), 
(right, bottom) ibidem, for = 8.8 TeV pPb (LHC). 

(right, top) The corresponding single inclusive pion spectra. 
The thin vertical line denotes the kinematic range with statis- 
tics of more than 1000 events per GeV-bin after one month 
of LHC operation with pPb. 

suggested (see e.g. [l3|]) that the characteristic enhance- 
ment of Rp A (pr,y) in the p*r-range of a few GeV (typi- 
cally referred to as Cronin effect [9j) can be understood 
in terms of transverse parton momentum broadening in- 
duced by multiple scattering. Transverse nuclear broad- 
ening is the prototype of a generic nuclear modification, 
for which we do not know whether and how it could be 
absorbed in collinear, process-independent nPDFs. How 
can one test whether the physics underlying Rp A (px,y) 
can be attributed to a nuclear modification of longitudi- 
nal parton momentum distributions and thus can indeed 
provide reliable quantitative constraints on nPDFs? To 
address this question, we have calculated Rppt,{PT, y) for 
the production of neutral pions in proton-lead collisions 
at the LHC, see the right hand side of Fig. Q] 

The LHC can collide protons and Pb ions with a max- 
imum center of mass energy of a/snn = 8.8 TeV. While 
pPb is not yet part of the initial LHC program, there 
are estimates 
nosity of C p pb 



14j that without major upgrades a lumi- 



10 



29 , 



could be achieved. With 



these assumptions, we find that running the LHC for one 
month would allow one to map out the single inclusive ir - 
spectrum up to transverse momenta well above px — 50 
GeV (see Fig.[T]). 

Remarkably, if the entire nuclear effect in pPb colli- 
sions can be factorized into nPDFs, then the shape of 
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the nuclear modification factor measured at the LHC will 
be qualitatively different from that observed at RHIC. 
This is so because at more than 40 times higher center 
of mass energy, final state hadrons at the same trans- 
verse momentum test O(40) times smaller momentum 
fractions Xi. As a consequence, Rp Pb {pT,y = 0) at the 
LHC will be dominated by the shadowing regime, and 
thus show a suppression, for pp < 10 4- 20 GeV, whereas 
RHIC data show a clear enhancement in this region. Fur- 
ther, LHC data will show a nuclear enhancement in the 
anti-shadowing dominated range of px ^ 10 4 20 GeV, 
whereas the RHIC nuclear modification factor starts be- 
ing dominated by ^-values in the EMC-regime. 

We emphasize that a shift of the maximum of 
Rp° Pb {pT,y = 0) to values of p T > 50 GeV at the LHC is 
a natural consequence of nuclear modifications in longi- 
tudinal parton momentum distributions, as encoded e.g. 
in EPS09. In contrast, no mechanism is known which 
could account for such a large py-shift in terms of trans- 
verse parton momentum broadening; the -^/i-dependcncc 
of transverse momentum broadening is much milder. The 
inverse is equally true: a mild shift of the maximum of 
R p°Pb(PT,y = 0)| v ^7 =8 .8TcV to values of p T < 10 GeV 
could not be accommodated naturally in a collinear fac- 
torized approach, since it would imply a nuclear enhance- 
ment of some PDFs below x ~ 0.01, which is inconsistent 
with the position of the anti-shadowing region. However, 
such a mild shift would be a natural consequence of trans- 
verse momentum broadening. 

We also emphasize that at the LHC, a collinearly fac- 
torized approach results typically in a mild enhancement 
of Rp Pb (j>T,y — 0) above unity for a wide transverse 
momentum range px > 10. In this kinematic range, sup- 
pression factors of order 2 are inconsistent with all ex- 
isting nPDFs. In contrast, models based on non-linear 
small-x evolution (see e.g. Ref. [15]) arrive naturally at 
such large suppression factors. 

The two examples mentioned above illustrate how the 
much wider kinematic range accessible in pPb collisions 
at the LHC allows one to discriminate decisively between 
qualitatively different models of nuclear modification. 

The dynamical explanation of the nuclear modification 
factor Rp Pb in terms of process- independent collinearly 
factorized nPDFs implies that, as a function of ^/snn and 
rapidity j/, the same nuclear effect manifests itself in very 
different kinematic ranges of pA collision data. As seen 
in Fig. [21 the rapidity dependence of Rp Pb allows one to 
scan the main qualitatively different ranges of standard 
nPDFs in an unprecedented way. At backward proton 
projectile rapidity (y = —4), where relatively large nu- 
clear momentum fractions x are required for hadron pro- 
duction, the nuclear effects in Fig. [2] are seen to be dom- 
inated by the anti-shadowing regime at low transverse 
momentum px < 20 GeV and by the EMC suppression 
at higher px- As one moves to larger rapidity, where 
smaller nuclear momentum fractions dominate hadron 
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FIG. 2: Rapidity dependence of Rppt, Q for ^/snn = 8.8 
TeV pPb (LHC). The different plots scan the dependence 
from y — —4 (close to Pb projectile rapidity) up to y = 4 
(close to proton projectile rapidity). Labels indicate whether 
the nuclear modification originates mainly from the shadow- 
ing (S), anti-shadowing (AS) or EMC regime. Vertical lines 
illustrate the rapidity-dependent px range, which can be ac- 
cessed experimentally with more than N events = 1000 (= 10) 
per GeV-bin within one month of running at nominal lumi- 
nosity. 

production, the anti-shadowing regime contributes up to 
increasingly high-p^, and opens up a wide window of 
transverse momentum, in which the shadowing region of 
nPDFs can be tested experimentally. 

Within the collinearly factorized approach, one expects 
non-perturbative corrections to the ansatz (TT]). How- 
ever, these corrections die out as inverse powers of the 
resolution scale. In contrast, while nuclear effects in 
nPDFs also depend on the resolution scale, their de- 
pendence is only logarithmic, so that sizeable nuclear 
effects are expected to persist at perturbatively large px- 
scales. Therefore, concise tests of the collinearly fac- 
torized approach require particle production processes 
at sufficiently large 'perturbative' momentum transfers. 
In pPb collisions at the LHC, the experimental access 
to a wide, nominally perturbative pr-range (jpx > 10 
GeV, say) is thus a qualitative advantage. In particular, 
the forward rapidity dependence of RHIC data @, [l6[ 
on Rj Au does not yet provide a decisive test for the 
collinearly factorized approach, since they test relatively 
low resolution scales, where large corrections to (P) could 
be expected even within the framework of a collinear fac- 
torized approach. For this reason these data have not 
been included in recent nPDF analyses [l|. In contrast, 
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FIG. 3: Rf Pb © for = 8.8 TeV pPb (LHC), from two 

different sets of nPDFs. 

within the framework of a collinearly factorized approach, 
one does not know of sizeable corrections to (U} in the 
range 20 < pt < 40 GeV, which will be uniquely accessi- 
ble at LHC and where characteristic rapidity-dependent 
features are seen in Fig. [5] 

So far, we have emphasized that well beyond quan- 
titative improvements, pPb collisions at the LHC have 
the potential to submit the very assumption of collinear 
factorization to decisive tests. In particular, a strong 
suppression of Rp Pb (pT,y = 0) at high pt > 10 GeV, 
or the persistence of the maximum of Rp P i,(pT,y = 0) 
at pt < 10 GeV is inconsistent with all current nPDFs 
and it tests an x-range for which existing data provide 
constraints. Therefore, if observed, such features would 
shed significant doubt on the use of the factorized ansatz 
((T|) for calculating nuclear effects, while they could be 
accounted for naturally in the context of qualitatively 
different dynamical explanations, mentioned above. 

Despite these perspectives for qualitative tests of 
collinear factorization, we caution that current global 
analyses of nPDFs come with significant uncertainties. 
While not all conceivable data on Rp Pb (pT,y) at the 
LHC can be accommodated within a collinearly factor- 
ized approach, a significant spread could. To illustrate 
this, we have compared in Fig. [3J the nuclear modifica- 
tion factor for two nPDF sets, which are known to show 
marked differences. In particular, in contrast to EPS09, 
the gluon distribution of HKN07 @ does not show an 
anti-shadowing peak but turns for x > 0.2 from suppres- 
sion to strong enhancement at initial scale Q 2 = 1 GeV 2 . 
Inspection of Fig. [3J reveals that for HKN07, the size 



and position of the maximum of Rp Pb (pT,y) at nega- 
tive y arises from an interplay between the nuclear en- 
hancement of the gluon PDF (which increases with x and 
hence with px) and the relative contribution of the gluon 
versus the quark distribution to Rh ph {pT,y) (which de- 
creases with pt). Fig. [3] thus illustrates that within the 
validity of a collinearly factorized approach, LHC data 
can resolve the qualitative differences between existing 
nPDF analyses and can improve significantly and within 
a nominally perturbative regime on our knowledge of nu- 
clear gluon distribution functions. Data on other single 
inclusive particle spectra and jets in pPb at the LHC can 
further constrain global nPDF analysis, thereby testing 
the concept of collinear factorization of nuclear effects 
and improving our knowledge of nPDFs as long as this 
test is passed. 
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